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A microscopic Hamiltonian theory of the FQHE, developed by Shankar and myself based on the 
fermionic Chern-Simons approach, has recently been quite successful in calculating gaps in Fractional 
Quantum Hall states, and in predicting approximate scaling relations between the gaps of different 
fractions. I now apply this formalism towards computing magnetoexciton dispersions (including 
spin- flip dispersions) in the u = ^, |, and | gapped fractions, and find approximate agreement 
with numerical results. I also analyse the evolution of these dispersions with increasing sample 
thickness, modelled by a potential soft at high momenta. New results are obtained for instabilities 
as a function of thickness for | and |, 



and it is shown that the spin-polarized | state, in contrast 



to the spin- polarized 3 state, cannot be described as a simple quantum ferromagnet. 
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I. INTRODUCTION 

The Fractional Quantum Hall (FQH) effectB has in- 
troduced us to new states of electrons in high mag- 
netic fields. Seminal theoretical progress was made by 
Laughlira, who showed that for filling fractions v = 
l/(2p + 1) the electrons form a strongly correlated state 
which is incompressible, and has quasiparticle excitations 
with fractional chargea e* = e/(2p+l). Subsequently the 
excitations were also shown to have fractional statisticsu. 

A unified understanding of all fractions v = pj (2sp + 
1) was achieved by the Composite Fermion picture of 
JairD. In this picture, the true quasiparticles are elec- 
trons dressed by 2s vortices, which are called Composite 
Fermions (CFs). At a mean field level, the CFs then see a 
reduced field B* = B/(2sp + 1), in which they fill p CF- 
Landau levels (CF-LLs), and exhibit the integer QHE. 
This picture has bejen very successful in obtaining excel- 
lent wave functionsEl. In the past few years, the experi- 
mental reality of CFs has also been firmly establishccB. 

Great progress has also been made in arriving at a 
functional integral description in which some of the non- 
trivial properties found in the wave-function approach 
arise at the mean field level. The functional treatment is 
based on the Chern-Simons(CS) transformation^], which 
performs fluxr-attachment viarthe CS gauge field to obtain 
either bosonsOtil or fermionstL3. These theories have pro- 
vided us with a link between the microscopic formulation 
of the problem and expjpr.irrient, both for incompressible 
and compressible statesHEJ. 

Recently R. Shankar and the present author developed 
a hamiltonian CS theory for the FQH-, statesESEJ. In- 
spired by the work of Bohm and PinesO on the 3D elec- 
tron gas, we enlarged the Hilbert space to introduce n 
high-energy magnetoplasmons degrees of freedom, (n also 
being the number of electrons) at the same time imposing 
an equal number of constraints on physical states. Upon 
ignoring the coupling between the oscillators and 
fermions we obtained some well known wavefunctionsH 



However the fermions still had the bare mass, and the fre- 
quency of the magnetoplasmons was incorrect. Hence a 
final canonical transformation was employed to decouple 
the fermions from the oscillators in the infrared limit. 

We choose to call the final fermions the composite 
fermions for the following reasons. Firstly, the final 
fermions have no dispersion in the absence of interac- 
tions and acquire an effective mass dependent on inter- 
actions alone. Next, the final canonical transformation 
assigns to each fermion the magnetic-^aoment e/2m as 
mandated by the arguments of RefsEjO. The central 
result of our formalism is the formula for the electronic 
charge density, which takes the following form, separable 
into high- and low-energy piecealj, at small q: 
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where A, A> refer to the annihilation and creation oper- 
ators of the magnetoplasmon oscillators, / = 1/y/eB is 
the magnetic length, and Hj = Pj + eA*(rj) is the ve- 
locity operator pf,the CFs. The oscillator piece saturates 
Kohn's theoremcd. The rest, to be called p, is obtained by 
adding to the canonicaiiy transformed electronic charge 
density a particular multiple of the constraintiLj (in the 
physical subspace, one can add any multiple of the con- 
straint without physical consequences, but we wish to 
work in the full space). It has some very useful proper- 
ties in the full space: 

• p satisfies the magnetic translation algebra 
(MTA)c3 to lowest leading order. Since this is the 
algebra of the electron density in the lowest Landau 
level (LLL), we have performed the LLL projection 
correctly in the infrared. 

• Note that p is a sum of a monopole with charge 
e* = e/(2p+ 1), which is the charge associated with 
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the CF, and a dipole piece which alone survives^ 
v = 1/2 and has the value proposed by ReadEI 
(A number of recent^onstructions have emphasized 
this dipolar aspectE3c3). 

• We also find that that as q — > all transition matrix 
elements of p from the HF ground state vanish at 
least as q 2 . 

The final property is an essential property of physi- 
cal charge density matrix elements from incompressible 
liquid ground states in the LLLc3. We will present argu- 
ments in the next section to show that if one intends to 
use the Hartree-Fock approximation ignoring constraints, 
these properties of p are essential. They make it plausible 
that p does not suffer vertex corrections. 

The Hamiltonian of the low-energy sector (dropping 
the magnetic moment term) is 



(2) 



where v(q) is the electron-electron interaction. We cir- 
cumvent the fact that p is to be trusted only for small 
q as follows. Consider real samples which have a finite 
thickness A of the same order as I, so that-the Coulomb 
interaction is cutoff at large w-ayevectorsEa. It was re- 
alized by Haldane and RezayiE3 that this has a large 
effect on the gap, while leaving the wavefunctions es- 
sentially unchanged. We will focus on such interactions, 
parametrized by A = A/ 1 for whichjiumerical results for 
the transport gaps are availableEatH The advantage is 
that as A becomes large only small-q matrix elements 
of the density are invoked in computing gaps, and we 
expect our theory to become more accurate. It is possi- 
ble that beyond some large A the incompressible liquid 
might cease to be the ground statecJ. Our theory, which 
is based on a liquid state with uniform density, can be 
expected to work up to this A. In fact, the magentoex- 
citon dispersions can be used to infer these instabilities, 
as will be shown later. 

This hamiltonian is to be supplemented by n con- 
straints which identify the physical subspace. We will 
expand on this issue in the next section. In two 
earlier papers, we presented calculations of gaps for a 
few fractional!! in the Hartree-Fock approximation, and 
tested certain scaling relations that arise naturally in the 
our theory against numerical results obtained using CF 
wave functions^. Emboldened by the good agreement we 
found between the predictions of our theory and numeri- 
cal results, in this paper I present results for magnetoex- 
citon (ME) dispersions for the polarized FQHE states at 
v = |, |, and j based on our formalism. I will also 
present results for the spin- flip MEs, one branch (not al- 
ways the lowest) of which corresponds to the spin waves 
in the ferromagnetic polarized states. 

The computation— of ME 

dispersions has a long-ijislorjO in the IQHE. Many dif- 
ferent approximationsLJta have been employed to obtain 
the ME dispersions in the IQHE, most of which can be 
subsumed into the powerful time-dependent HF (TDHF) 



treatment of MacDonaldcJ'tS. In this approach one al- 
lows the hamiltonian to scatter particle-hole excitations 
with different Landau-level indices into each other, and 
diagonalizes the resulting matrix. 

A serious problem arises when considering ME disper- 
sions in the FQHE. Since the states are strongly corre- 
lated in terms of electron variables, it is difficult to find 
good approximation techniques. Consequently, the most 
trusted results-are those of exact diagonalizations on fi- 
nite systemsEHEZI. There are also results from other methr- 
ods, such as the single-mode approximation (SMA)E3, 
which works well near the roton-like minimum for v =rk x , 
and the use of CF wavefunctions for excited statesEa. 
Field-theoretic approaches include the seminal work of 
Lopez and Fradkinta for the collective modes of gapped 
fractions based on a fermionic CS approach, and an RPA 
treatment^ based on the formalism of Halperin, Lee, and 
Read (HLR)tLj. These field-theoretic approaches suffer 
from the problem that the bare mass (actually the band 



mass m&)E3 or a phemenological mass m* enters the dis- 
persions, whereas in the LLL the scale should be set solely 
by e 2 /el. 

Our approach essentially solves the problem of strong 
correlation by rewriting the theory in terms of CFs, which 
are the quasiparticles of the theory in the same sense as 
those in Landau's Fermi Liquid theory: The quasiparti- 
cles can interact strongly with a Fermi-Liquid-like inter- 
action, but matrix elements which scatter them out of 
their states are small at low energies, allowing them to 
be long-lived. Furthermore, the only scale in our low- 
energy hamiltonian is indeed e 2 /el, as desired. For the 
reasons to be explained below, we can use a formalism 
such as TDHF, which ignores vertex corrections, with our 
formula for p. 

It must bej—Mnphasized that the TDHF approach 
of MacDonalcOEil which includes many Landau levels 
(LLs) is absolutely essential for treating the FQHE case. 
In the integer QHE one has two distinct energy scales. 
Tilo c and e 2 /el, and one may go to the strong-field limito 
fiLo c >> e 2 /el to make the problem theoretically conve- 
nient without losing any essential physics. In this limit 
only the lowest unoccupied and highest occupied LL need 
be considered, and LL- mixing of the MEs can be ignored. 
However, in our LLL analysis of the FQHE, there is only 
one energy scale, and the variation of energy within an 
ME mode is comparable to the splitting between differ- 
ent modes. We will see that the correct magnetoroton 
minimum cannot be obtained without considering this 
"CF-LL mixing" . 

The spin quantum number produces additional non- 
trivial features. It was realized by Sondhi et ac3 that 
even in the absence of Zeeman interactions, the v = 1 
state, and correspondingly the Laughlin fractions l/2p + 
1, were ferromagnetic due to exchange interactions. The 
symmetry-breaking leads to gapless spin waves with a 
quadratic dispersion. If the Zeeman interaction is now 
turned on the spin waves acquire a gap. Due to the fact 
that all other excitations of the system are much higher in 
energy than the spin waves, the v = 1 and the Laughlin 
fractions are dynamically identical to a simple quantum 
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ferromagnet at low temperatures^. The detailed predic- 
tions have-been dramatically confirmed by experimen- 
tal resultsc3. There have been few attempts to compute 
the low-energy spin-flip modes for the polarized FQH 
states. One, due to Nakajima and Aokia, relates the 
FQH spin-wave mode to the IQH spin-wave mode, per- 
forming the flux-attachment by shifting the pseudopotea.- 
tials. This, when tested against exact diagonalizationto, 
produces good results for the v = ^.^and \ spin- wave 
modes. Another approach, by Mandallil, uses fermionic 
CS theory to compute the response functions, and from 
them the spin-flin-dispersions. In the spirit of the earlier 
CS approachesLL3'E3 a phenomenological effective mass is 
used to approximate the energies of the CF-LLs, and the 
strong-field approximation is used. Where appropriate, 
we will comment on the relationship between these and 
our results. 

The plan of the paper is as follows. In section II we 
will set up the conventions and define various terms. We 
will also show that the states with p filled CF-LLs are 
indeed HF ground states for v = p/2p + 1, and we will 
make some important remarks about the applicability of 
TDHF in our approximation scheme for the gapped frac- 
tions. In section III we describe the method of comput- 
ing the ME dispersions, which is essentially MacDonald's 
approaches expressed in an operator form. In Section IV 
we present the results for the spin-polarized and the spin- 
flip ME dispersions and compare them to previous work. 
As mentioned earlier, beyond some A the incompressible 
FQHE state becomes unstableEJ. Other instabilities ia- 
clude the transition to a Wigner crystal at small fillingO. 
The latter instability has tepa studied by looking, at the 
ME dispersion in the 1/7l3u3 and 1/9 systems^. It is 
found that -the ME dispersion becomes negative at finite 
wavevectoiO for the 1/9 FQHE state, indicating that the 
ground state is unstable to the creation of large numbers 
of MEs. In analogy with this work, we will go back to the 
■| , | , and | FQHE states at large A and find instabilities 
in the ME dispersions for v = | and |. We end with 
conclusions and open questions in section V. 



II. SETTING UP THE FORMALISM 

Let us return to the hamiltonian of equation (2), fo- 
cussing for simplicity on spin-polarized states. Since 
we have a many-body Hamiltonian, we should choose a 
many-body ground state (GS). The formula for p con- 
tains the velocity operator for the CFs in the effective 
field B* = B/2p+l, suggesting that the CFs fill CF- 
LLs. 

We now have a problem where the states are the 
same as those for (spinless) fermions in the IQHE. We 
work in the symmetric gauge, in which the single-particle 
wavefunctions are characterized by the LL index n and 
the angular momentum index k. The-formalism in this 
gauge has been extensively developedacZI in the context 
oLanyon superconductivity. To conform to the notations 
oicZl we will define the magnetic field to be along the pos- 



itive z direction, which produces the wave functions 
V2TT2 n+k nlk\ 2 2 

(3) 

Here z = x + iy, and all distances are measured in units 
of the effective cyclotron length I* — l\J2p + 1. Note that 
the LLL wave function is antianalytic with this choice of 
field direction. 

The density is a one-body operator, and can._ereate 
particle- hole pairs above the GS. It is well-knownc3 that 
these pairs come with a conserved momentum q, which 
also characterizes the separation between the particle 
and the hale r p h = z x ql* 2 . The magnetoexciton wave 
functionsEj ip nn > (q; r, r") describe a particle in the r » t ' 1 LL 
and a hole in the n nh LL. In the notation of refj£3, it is 
expressed as 

vw(?»") = ( -^-, — =r—S 2d i - f ) n (2& - i )"' 

V2tt2" + " n\n'\ 

e -i(ri+ri+Q 2 ) e ^(z* 1 z 2 +z* 1 z q ~z*z 2 ) 

Here L is the linear size of the system in units of I* (the 
dimensionless area is L 2 ), and z q = iQ+ = il*{q x + iq y ) 
is the dimensionless wavevector. 

We can now define the ^-coefficients, which denote the 
overlap between particular single-particle + single-hole 
states and the ME wavefunction as (for k > k') 

= (-iT'^^{^) k ~ k 'L k k r k '(y)e-y/ 2 (5) 

Here is a Laguerre polynomial, and y — Q 2 /2. The 
x-coefficients are also the matrix elements of the uni- 
tary transformation that relates the product particle- 
hole basis with the ME basis. In order to obtain the 
x-coefficients for k < k! we use 

x k :, k n (q) = (-lT +n '(x k £(-q)y (6) 

Now we are ready to express p in this basis set. The 
form of p presented in Eq.(|]) is correct for the field di- 
rection being —z. In order to conform to the notation we 
are usingc3 we have to use the following form 

p(q) = ^+7 + ^ 2 £> x i,,X3 ) ( ? ) 

j 

where H,- is now the velocity operator in the effective field 
pointing in the +z direction. Now we invert Eq. (^|) to 
get 



d (^1)^2(^2) = 




(9))*V'n 2 n 1 (g;r2,ri) (8) 



where we have used the compact notation vi = (m, ki). 
To express the density in terms of second quantized op- 
erators we note that 
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We begin by writing down the hamiltonian in terms of 
the CF creation and annihilation operators 



5j Pn 2 n 1 {<l){x l i 2Vl 



(9) 



where we define 



(16) 



/Wi (<z) = J d 2 rt/j n2ni (q; r, r) + H 2 q x n e 



l h^T ( -iQ+ \ n 2-m 

2p+l 73^ V "2^ \/2 / 



Consider the transition matrix element of H be- 
tween two different single-particle states p = (to, k) ^ 
(to', fc') = p' on top of the assumed GS of p filled CF- 
LLs; 



xe-y/ 2 (n 2 L n n lzr + 2L%-*» - {m + l)L^) (10) 

This last expression is valid for > n\. For n-i < ni 
one uses the relation 



(11) 



In what follows, it will be useful to define the magnetoex- 
citon operators 



Xmmio) = zlix^iq))* dl^ 



(12) 



k\k 2 



where the sums are over the angular momentum indices. 
Finally, one can compactly express the electron density 
operator as 



(0) 



(13) 



Apart from the trivial dependence of x on n' , we see 
the separation between the angular k labels and the "ra- 
dial" n labels. 

We will need the following identities, which can be 
easily established by using the defining relation for x, 
Eq. (|); 

= (-1)" ^ % (I 4 ) 

k 

Y,<^ X ™M) = (-l)»^eT<3i*$ a x 
k 

xtttii+fc) (I 5 ) 



A. Calculation of Gaps 

As mentioned before, we have to assume a many-body 
ground state. Since the low-energy hamiltonian contains 
the velocity operator in the effective field B* = B/2p+l, 
it is natural to assume a ground state of p filled CF-LLs 
for v = p/2p + 1. This is not an eigenstate of H, which 
can create particle-hole pairs above it. However, it is the 
HF ground state. We now establish that, for rotationally 
invariant interactions, the Hamiltonian does not mix a 
single-particle (or single- hole) state with any other single- 
particle (or single-hole) state. This is the signature of a 
HF state. 



< GS\d^Hdl,\GS >=\J ^L v (q)J2 

Pntm (-q)Pn 3 n 2 (tf)xt iVl (~Q)K 3 u 2 (?) 

< GS^dl^A^dl, \GS > (17) 



Now we perform the usual Wick contractions (at equal 
time) with 

< GS\4d v '\GS >= 5 vu ,N F {v) (18) 
< GS\d u dl \GS>= <W(1 - N F (v)) (19) 

where Np{v) denotes the occupation of the single- 
particle state v. In performing the Wick contractions 
one should take care not to contract the two operators 
belonging to the same density (say d\, with d„ 4 ), since 
this leads to the q = density, which is cancelled by 
the background positive charge. Now we can write the 
transition matrix element as 

(1 - N F (p))(l - AW))/ (Cf^^E ( 
(1 - N F{v))xl tl {-q)x* l , lJ {q)p nm {-q) Pmln {q) - 

N F(^)xl^(~q^xl^p m , n (~q)p nm (q^j (20) 

We now use the x-identities to perform the sum over 
the angular momentum index of v in the two terms to get 
a factor of x* , (q = 0) in both terms. The explicit form 
of x tells us that setting q = forces k — k' . This leaves 
the integral over the reduced density matrix elements. 
Consider the first term for specificity 



tt nS kk> f d 2 q \ - 

Hl = —J (27T (9) ^ ( ~ 1 

n>p 



Pmn(-Q)Pnm'(q) 



(21) 



By using the explicit functional form of the density 
matrix elements, Eqs. (10) and (|Ti"|), we can reduce this 
to the form 



7T(5fefc' f d q ^ , 2 



jj&im—m 1 ) 



(22) 



where 6 is the angle of q with respect to the q x axis. It 
is then clear that for to 7^ to', the rotational invariance 
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of v(q) will force this term to be zero. The same result 
can be verified for the second term, and also for the case 
of two different single-hole states. 

This establishes the HF nature of the GS, and supports 
our view that we are working with the right variables. We 
can now verify from Eq. ( |10[) , as claimed in the introduc- 
tion, that as q — ► all transition matrix elements from 
the HF ground state vanish at least as q 2 , an essential 
property of charge density matrix elements in the LLIxj. 

Let us now find the energy of a single-particle state 
above the filled sea. This is just the diagonal matrix 
element of H, and we have already done most of the 
work. Setting m = mf in the above, we get 

< m ) = ^J ^v{ q )^{l-2N F (n))\ Pnm {q)\ 2 (23) 

The first term is the Hartree energy while the second 
represents the exchange, or Fock part. The hamiltonian 
can be normal ordered to make this obvious in operator 
form 

H= H Q + V 

Ho= £eo(n)dU< (24) 

V 

X Pn 4 n 1 (-9)Pn 3 n2(9)4i4 2 d V3 d Ui (25) 

eo(n)= M ^v(q)j:\Pnm(q)\ 2 (26) 

n 

Similarly, one can easily derive the hole energy, which 
is identical to the above, except for a minus sign. Notice 
that the energies depend only on the CF-LL index, and 
not on the angular momentum, and that there is a sum 
over an infinite tower of CF-LLs in the energy. 

So far our discussion is applicable to arbitrary v(q). 
Below we will find it convenient to specialize to a model 
potential that approximates the effects of sample thick- 
ness v(r) = e 2 / 'eV r 2 + A 2 , whose Fourier transform is 
v(q) = 2ne 2 exp (— qA)/eq. The parameter A is related 
to the sample thickness. 

At this point, two important questions arise: We have 
treated the hamiltonian in the HF approximation, and 
we have ignored constraints altogether. It is important 
to understand whether these approximations are reason- 
able, and how they might (at least in principle) be im- 
proved. The next four subsections will address these is- 
sues in turn. 



B. Why is Hartree-Fock reasonable? 

Hartree- Fock is an approach that works best when cor- 
relations are weak. The ground state is sought to be ex- 
pressed in terms of a single determinant of single-particle 
wave functions. Clearly, carrying out a HF calculation 
in terms of electron coordinates would be nonsensical in 
the EQHE problem: This is the message that Laughlin's 
worm has driven home. However, what we are doing here 



is radically different. We are working in terms of CF co- 
ordinates. The correlated ground state of electrons can 
be viewed as the ground statejpf independent CFs, which 
is the message of Jain's workcl. In other words, the sys- 
tem of electrons reorganizes itself into a state where the 
CFs are the Fermi-liquid- like quasiparticles. Thus one 
can reasonably expect to get good results from a simple 
approximation such as HF when working in CF coordi- 
nates. 

However, one other condition has to be satisfied if the 
results are to be numerically accurate. We have to make 
sure that vertex corrections are small. Large vertex cor- 
rections imply that the particle we started out with is 
not the physically observed quasiparticle, i.e., that the 
observed charge, dipole moment, etc are different from 
those entering the microscopic hamiltonian. In this case, 
the HF excited states will have small overlap with the 
true excited states, and the HF energies will be corre- 
spondingly poor. Our formula for p does have the correct 
charge and dipole moments of the final physical particles. 
In other words, they are already fully dressed, and ther- 
fore vertex corrections are expected to be small. 

There is another way to conclude that vertex correc- 
tions are small for p: RecallliS that in order to get p we 
took a particular multiple of the constraint and added 
it to the "naive" transformed electron density. In prin- 
ciple, one can add any multiple of the constraint to the 
density, and the physical matrix elements should be unaf- 
fected. However, this particular combination has a very 
nice property which we have mentioned before: All its 
transition matrix elements from the HF ground state are 
order q 2 or higher. This is a necessary property of all 
physical .density matrix elements in the incompressible 
fractions^. This property does not hold for the naive 
transformed density, or for any other linear combination 
of the density and the constraint; they all have transition 
matrix elements of order q from the assumed HF ground 
state. Of course, if the constraint were exactly imple- 
mented, it would force all physical matrix elements to 
behave correctly by producing vertex corrections. Thus, 
every combination except our p is guaranteed to suffer 
strong vertex corrections. Therefore we come to an im- 
portant conclusion: If we intend to use HF without im- 
plementing the constraint exactly, we must use p, and 
only p, for the electron density operator. 

C. Constraints 

Since it is very hard to implement these constraints 
exactly, we will use two different approaches. The first 
approach is to ignore the constraints altogether. The ra- 
tionale for doing this is the following: Suppose we had the 
exact expression for the electron density operator in final 
coordinates. Since this operator is gauge invariant {i.e. 
commutes with the constraints), it will have no matrix 
elements between the true physical ground state and un- 
physical states. Thus the sum over physical states can be 
extended to a sum over all states with no problem. How- 
ever, this requires that we also have the correct physical 
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ground state. In practice we have neither an exact ex- 
pression for p e nor the true ground state. However, our p 
is close to an operator that is gauge invariantea (we will 
elucidate this in subsection E) , and we can hope that the 
ground state we choose is close to the true one, so that 
all the above statements are approximately correct. 

While this approach indeed produces reasonable an- 
swers for the thickness parameter A > 2.0 or so, it pro- 
ducesnap unphysical divergence of the transport gap for 
A = fB. The reason is not far to seek. For small A the 
potential has significant weight at large q, where our ex- 
pression for p is not to be trusted. Thus in this region, 
our p couples to unphysical states, of which there are an 
infinite number, which produces the divergence. Thus, 
in order to get a physically sensible answer for small A 
one is forced to account for the constraint in some way. 
We choose the simplest possible way, by cutting off the 
high energy CF-LL states. There are many ways of see- 
ing that such a cutoff makes physical sense. Firstly, for 
a finite number of particles, there are a finite number of 
states in the many-body Hubert space (after projecting 
to the LLL, which we do by freezing oscillators), equiv- 
alent in dimensionality to the Hilbert space generated 
by 2p + 1 CF-LLs. Secondly, density matrix elements 
in the LLL are suppressed at high momenta by gaussian 
factors. Matrix-elements of our p between higher CF- 
LL states are peaked at higher momenta, and thus it is 
plausible that the higher CF-LLs are less physical. Fi- 
nally, even in the LLL, we expect CFs to describe only 
the low-energy dynamics. At .energies of the order of the 
first Haldane pseudopotentiaEJ V\ , we expect the CF to 
break up into electrons and vortices. Thus, CF-LLs can- 
not be the correct description at high energies. On the 
other hand, there is ample numerical evidence that the 
low-lying CF-LLs are very close to true physical statesEil. 

Thus we have the following picture: The higher CF- 
LLs probably have very little overlap with physical states, 
while the lower ones have large overlap. This argument 
does not tell us precisely which states we must keep. In 
the rest of this paper I will choose a sharp cutoff in the 
CF-LLs, wherein some number (around 2p + 1) of the 
lowest ones are kept, and the rest discarded. We will 
also investigate how sensitive our results are to a change 
in this cutoff, and we will find that for A > 2.0 or so, 
the results become insensitive, and that one can keep 
all the CF-LLs without much altering the results. The 
main operational aim of this method of implementing the 
constraints is to obtain qualitatively reasonable (though 
perhaps numerically untrustworthy) results for small A. 

Thus far we have given arguments mostly from within 
our theory. Of course, the correct way to settle mat- 
ters of principle is to implement the constraints from the 
start on an equal footing with the hamiltonian. Such a 
calculation has pot been carried out for the gapped frac- 
tions, but ReadE3 has recently studied the case of v = 1 
bosons in detail (this belongs to class of compressible 
systems, panaang which is included the half- filled LLL for 
electronaij't-j) , where constraints-are CT 1 ] r ' ig1 in restoring 
compressibility to the systemOOaOEij. We devote the 
next subsection to the lessons we can draw from Read's 



wor 




D. Read's Conserving Approximation for v = 1 
Bosons 

ReadH starts from the LLL formulation of Pasquier 
and HaldaneEa for v — 1 bosons, and implements 
the constraints, (which, unlike the constraints in our 
formulationE_3'c3, do not commute among themselves, 
but form an -algebra) in a path-integral conserving 
approximatiorO. Our greatest interest is in the gaps 
and ME dispersions, which can be read off as the poles 
of the density-density correlator. Let us therefore focus 
on his results for the irreducible density correlator Xirr, 
for which he findso that: 

• Xirr breaks up naturally into two additive contribu- 
tions. The first is the density correlator of dipoles, 
while the second represents the dipoles interacting 
with a transverse gauge field. It should be empha- 
sized that the fermions that enter the lagrangian 
have charge e, and that a nontrivial renormaliza- 
tion has happened due to vertex corrections in or- 
der to make them dipolar. 

• The contribution due to the transverse gauge field 
dominates at very low frequency, restoring com- 
pressibility to the system. This contribution con- 
trols the to and the w" 1 moments of the density 
spectral function. 

• At higher frequencies, the purely dipolar part dom- 
inates. For example, the uj moment of the density 
spectral function is identical to the result for purely 
dipolar fermions. 

Let us view our method in the light of these results. Our 
fcrmionic density, written in terms of p, is dipolar (for 
v = j) already at tree level, and no further vertex cor- 
rections are necessary. What is missing in our approach 
is the coupling to the low-energy transverse gauge field. 
However, this contribution seems to be relatively unim- 
portant at high frequencies. 

While we are not certain how all this might extend to 
the gapped fractions, we can make the following plausi- 
ble guesses: The vertex corrections will renormalize the 
charge of the fermion from e to e* , and produce a dipole 
moment (and perhaps higher multipole moments). Once 
again this renormalization (upto dipole terms) is present 
in our p already at tree level. Secondly, we can expect a 
gauge field to enter the picture. However, we once again 
expect it to be important only at low frequencies. Since 
all the answers we seek in the gapped fractions are at 
finite frequencies of the order of the gaps, we do not ex- 
pect the gauge field to make a large contribution to our 
results. 

To summarize, the effects of the constraint in Read's 
calculation are twofold: They renormalize the charge and 
dipole moment, and they create a transverse gauge field 
to which the dipolar fermions are coupled, and which is 
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important at low frequencies. The first effect is present in 
full in our approach, while we argue that for the physical 
quantities of interest to us, the contribution of the second 
can be safely ignored in the gapped fractions. 

We now turn to a very recent new formulation for the 
gapped fractions. 



E. Shankar's New Formulation for Principal 
Fractions 

Very recently, ShankarE] has proposed formulas for the 
charge density and constraint for all principal fractions. 
He writes both the electronic charge density (call it ps) 
and the constraint (call it xs) as a sum over exponentials 
of single-CF coordinates. The first two terms of each ex- 
ponential are exactly the two terms derived in our earlier 
workEj. ps and xs have some very nice properties: 

• ps satisfies the MTA to all orders in q. This for- 
mulation is not restricted to the infrared. 



state. Shankar points out@ that if one is interested in 
computing numbers in the gapped fractions, this second 
formulation, with many nontrivial renormalizations built 
in at tree level, may produce better numbers than the 
conserving approach. 

This second viewpoint meshes perfectly with our phi- 
losophy. It also gives us an algebraically consistent way 
to continue the expression for p to all orders in q. When 
one looks at the corrections to our p by comparing it to 
Shankar's p pre f, one finds that the leading order correc- 
tion (order q 2 ) vanishes, and the next|-cprrection (order 
q 3 ) is down by a large factor (w 500)ej. This suggests 
that while our p was derived at small q, it regime of va- 
lidity may be larger than we have any right to expect. 

To summarize the past four subsections, we have pre- 
sented several lines of reasoning that suggest that by us- 
ing our p for the electronic charge density, we are work- 
ing in a preferred representation where the effect of con- 
straints is minimal, and approximations such as HF and 
TDHF can be expected to give good answers. We now 
proceed to the calculations. 



• xs also satisfies the MTA with a modified mag- 
netic length, and is identified as the vortex charge 
densitycj. (It should be noted that this is different 
from our prigjnal formulation, since our constraints 
commutetj'E3, and should continue to do so after 
the canonical transformations). 

• ps and xs commute with each other, so ps is gauge 
invariant. 

• The hamiltonian can be written in terms of ps, 
which makes it gauge invariant, and amenable to a 
conserving approximation. 

Such a conserving approximation would be the way to 
settle matters of principle, such as the quantum num- 
bers of the excitations in the gapped fractioBSi and their 
statistics! However, Shankar also points outc2lthat there 
is another possible formulation of the problem that is 
very closely related to that followed in this paper, and 
provides a justification of it. He notes that one can form 
a preferred charge density 



III. METHOD OF CALCULATION 

Conventionally, in using the TDHF approximation^!, 
one first computes the density matrix under a time- 
dependent perturbation, and thence the behavior of the 
density with time to linear order in the external potential. 
The susceptibility is directly connected to the density- 
density correlator, and its poles give the ME dispersions. 

We will be following an operator variant of this 
method, where we will directly determine the time de- 
pendence of the ME operators (defined in Eq. (|l2|)) by 
commuting them with the hamiltonian. This method also 
has an ancient historjCJ. The expression on the right 
hand side will, in the spirit of HF, be simplified by re- 
ducing four-fermi operators to two-fermi ones by taking 
averages. 

We will initially consider the simpler spin-polarized 
case, and later show how to extend the formalism to in- 
clude spin. 



Ppref = PS + CXS 



(27) 



where c is the charge of the vortex. This charge den- 
sity has a very appealing physical interpretation, since 
each composite fermion has an electron and a vortex in 
close proximity, and p pre f is just the total charge of such 
a composite. One can now express the hamiltonian in 
terms of p pre f, and it commutes weakly with the con- 
straints ([-H,x] — x)- In the physical subspace this is as 
good a formulation as the one in the previous paragraph. 

Most important for our purposes, the first two terms of 
this preferred charge density are exactly the two terms in 
our p, Eq. (Q). Therefore p pre f enjoys all the nice proper- 
ties of p, such as explicitly displaying the right charge and 
dipole moment, satisfying the MTA at lowest leading or- 
der, and having q 2 matrix elements out of the HF ground 



A. Spin- Polarized Magnetoexcitons 



For a Heisenberg operator we have the relation 



(28) 



Let us begin by considering the time dependence of the 
ME operator X mi „ l2 (q) by computing [H, X mim2 (q)}. 
We will find it convenient to use the form of H split into 
Hartree and interaction terms as given in Eqs. ( p4|p5| ). 
It is then easy to verify that 

mim2 

(q) (29) 
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where eo is just the Hartree energy of Eq. (E 
compute the next term we will need 

[d Ul d V2 d Vz d L , i , djjj d^ 2 ] 



&n iv i d\j l d\, 2 d U3 d^ 2 
~ ^ ^iu 3 d L , i d l , 2 d L i 4 dfj l2 



3fj.2»i difii d\ 2 du 3 d Vl 
^2 "2 d Vi d Vz d Vl 



% To 



(30) 
(31) 



The hamiltonian has the symmetry (yi, ^3) — > (^2, ^4), 
and hence in the following we will keep only the first two 
terms (Eq. (|3(])), while removing the factor of \ in front 
of the interaction. Also, in TDHF one has to reduce 
all the four-point to two-point terms. We do this by 
contracting one creation and one annihilation operator in 
all possible ways, remembering that the operators that 
belong to the same density are not allowed to contract 
together (such a term would correspond to a uniform 
density, which is cancelled by the positive background). 
These terms can be grouped as follows: 

N F {mx)8 IJ , 2Vl {S IJ , lV3 dl 2 d Ui - Sfj, lUi dt 2 d V3 ) (32) 
v\ti?.d V2 d Vz 8 V2jl2 dt^ (33) 
+ N F (n 2 )6 l _ l2 ,y 1 5v 2 „ i dl li d„ 3 - Npin^S^^S^d^d^ (34) 

Now we multiply the above by the appropriate factors 
and integrate over <f , and use the x-identities. It is easy 
to show that the last line (Eq. ([34])) produces just the 
Fock energies of mi, TO2, which occur in the right combi- 
nation to add to [Hq,X] to give 



(e(mi) - e{m 2 ))X mim . 2 (q) 



(35) 



Let us now consider the more interesting terms from 
Eq. ( |32] ) and Eq. (|33"|). There are four such terms, and 
they are equal in pairs. Performing the Wick contractions 
and using the identities (bearing in mind that 6^ — 



(2*r 
~1T 



S 2 (Q) 



_ (27T) 



L'i\i*yi "{<l))i we obtain 



[V,X mim2 (q)} = - £ X niri2 {q)(N F ( mi ) - N F (m 2 )) x 

mn.2 

((-l)" ll+m2 W ^v(q)p mim2 (-<f) Pn2ni (q) 



(38) 



The hamiltonian can now be thought of as a matrix in the 
space of "naive" MEs, which are represented by the X 
operators. These are not eigenoperators of H, which has 
transition matrix elements between different naive MEs. 
An important point is that both positive and negative en- 
ergy naive MEs are connected by the hamiltonian, which 
is not hcrmitian. Nevertheless, the eigenvalues will occur 
in pairs ±Wj, with ujj > being identified as the energy 
of a ME mode. Finally, notice that both terms of an ar- 
bitrary transition matrix element go to zero in the limit 
of large q, the first because the matrix elements p nn < are 
suppressed with gaussian factors, and the second due to 
the oscillatory factor e l Q xS _ Thus the naive MEs decou- 
ple and become the eigenstates of the TDHF hamiltonian 
in the limit of large q. 

The density matrix elements p ni n 2 are complex, which 
means that generic Hamiltonian matrix elements are also 
complex. To simplify the future discussion we define the 
following reduced density matrix element p by the defining 
relation 



Pnyn 2 {q) 



L 



1 



(2p + 1) 



e i( ni -n 2 )(e q -n /2)p nin2 (Q) (39) 



where p is real, symmetric in nin 2 , independent of the di- 
rection of q, and is a function of Q — I* q alone. It is clear 
that the p matrix element can be made real by choosing 
9 q = tt/2, thus making the direct matrix element of H, 
Eq. (pq), real. The same choice also makes the exchange 
matrix element, Eq. (|37|), real. The fact that all matrix 
elements of H in the ME subspace are real speeds the 
numerical diagonalization. 



B. Scaling of ME Dispersions 



(36) 



where Q = l*q, S = l*s. For the case of the diago- 
nal matrix elements, i.e., {n\n 2 ) — (mim 2 ), these two 
terms have a simple physical meaning. The first term is 
the direct interaction between the particles, which shows 
up as an exchange interaction between the particle and 
the hole, while the second is the exchange interaction 
between particles, or the direct Coulomb interaction be- 
tween the particle and hole. We can now write the entire 
commutator as 

[H,X mim2 (q)] = (e(mi) - e(m 2 ))X mim2 (q) 

- (N F ( mi ) - N F (m 2 )) £ (-l) mi+roa X nina ($) x 
nin 2 



The discussion so far is valid for an arbitrary potential. 
Let us now specialize to our model potential 

v(q) = e~ Aq = e" A Q (40) 



Q 



where A* = hjl*. Using the definition of the Bessel func- 
tion 



Jn{Z)- j 



(41) 



and fixing 8 q — 7r/2, we can perform the angular in- 
tegration and cast the direct and exchange matrix ele- 
ments into the form (ignoring the factor of —(N F (mi) — 
N F (m 2 ))) 
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V d ( mi m 2 : nina ) = if p ^ffl«|!iffiM e A*g ( 42 ) 
^(mmu : mru) = -A^(-l)™ 1+ ™ 2 x 

OO 

(QS) (43) 



Here A%, 



i p — ; ( 2 p+i)5/2 ^ S a common energy scale. It turns 
it the fac 
energies to get 



out that the factor K p can also be extracted from the HF 



e(m) 



e(m) 



OO 

\ I dSj2\Pnm(S)\ 2 (l~2N F (n)) (44) 



Since the overall factor of K p can be extracted from 
the entire TDHF hamiltonian, it follows that the ME dis- 
persions scale with K p . Now suppose we had analyzed 
the fraction v(jp,r) — p/(2pr + 1). Everything above 
would have followed unchanged except for the replace- 
ment 2p+l—>2pr+l, and hence K p and A* would change, 
but nothing else (it is important to note that the ground 
state still consists of p filled CF-LLs). Thus, in our ap- 
proach, the ME dispersions at the same p, but different 
r, are related by scaling. This is an extention of the 
scaling relatipns for gaps that we verified in an earlier 
publication!! 2 ]. This scaling holds exactly in our theory 
provided the constraints are ignored. If the constraints 
are taken into account by cutting off the CF-LLs, then of 
course a r dependence will creep into the ME dispersions. 



C. Spectral Function and Spectral Weight 



It is frequently of interest to determine how much of 
the spectral weight is being carried by various modes, 
and in particular, the lowest mode. Such considerations 
are crucial for the success of the single mode approxima- 
tion (SMA)t3, which is known to be good near the roton 
minimum for |. Let us go back to the TDHF hamilto- 
nian of Eq. (|3^) , and assume that we have diagonalized it 
into modes a 



ttTDHF (fAxT/Ra /-a 
n\U2 W/ nin2:mim2 vi) 



(45) 
(46) 



where the sum over nin2 is implicit. Note that since 
jjTDHF j g n0 £ nernu tian we need to consider both right 
and left eigenvectors. 

Now we can go back and add an external time- 
dependent potential V (q, uj)e~ lult that couples to the elec- 
tron density in the hamiltonian. It is easy to verify that 
the equation for the response of the magentoexciton op- 
erators to the external potential is 



(u>5 rnim2:nin2 + H mim2 . nin2 )X nin2 (</, ui) — 
{N F {m x )-N F {m 2 )) ^£/Wh(£) (47) 

We can resolve any function of H TDHF [ n terms of its 
eigenvectors and eigenvalues, and in particular, we can 
invert the above equation to obtain 



a:niri2 

(N F (m) - N F (n 2 )) 



l* 2 L 2 



u) + E a 

(q)V(q,u) (48) 



Recalling the expansion of the density in terms of the 
magnetoexciton operators, Eq.(|l3|), we obtain for the 
density response 

(P)= E Pm.mM^mM^Er<%M 
* — * " UJ -+- £/ a 



a:mim27T-in2 



(N F {m) - JV>(n 2 )) , » t a / * \ 

x ~ i«2 L 2 - pn 2 m(q)V(q,uj) (49) 

Now we use the fact that due to the peculiar structure 
of H TDHF for spin-polarized MEs *^% 2 = {N F {mx) - 
^('W2)) , l l m" m2 - Further, using the definition of the re- 
duced density matrix elements, Eq. ([59]), we obtain the 
susceptibility to be 



A a 



uj + E n 



(50) 



where 



E Pm 2 mM^mM 
mirn-2 

2n(2p+ l) 2 /* 2 



(51) 



In the TDHF approximation all the modes are sharp, 
and do not decay. The spectral function is then a se- 
quence of delta functions at the energies ±E a with weight 
A a . Of course, if we go beyond TDHF, all the modes 
except the lowest will acquire a finite width. Even the 
lowest mode may have a dispersion which kinematically 
permits a quantum of high q to decay into two or more 
quanta of lower q. 

To summarize, our approach to computing the spin- 
polarized ME dispersions will be to compute the TDHF 
hamiltonian matrix in the subspace of MEs, and to di- 
agonalize it. If constraints are ignored, this subspace is 
infinite-dimensional, since there is an infinite tower of 
CF-LLs. We will follow the usual prescription of increas- 
ing the number of CF-LLs kept till the answer converges. 
For small A we find that the gap is unphysically large, 
and diverges as A^O. Here we will follow the discussion 
presented earlier and keep a finite number of CF-LLs, 
approximately equal to 2p + 1, to get a qualitative idea 
of the predictions of our theory in this regime. It should 
be emphasized that we trust the results only for A > 1.5 
or so. Fortunately, this seems to-be-close to the physical 
regime for high-density samplesEiLlj. 

Let us now turn to the spin- flip MEs. 



D. Spin-Flip Magnetoexcitons 



In order to incorporate the spin quantum number we 
must go back to the beginning and ask how the flux at- 
tachment is accomplished in the initial representation. 
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As implicit in the early wave-function formulation of 
HalperincSI, and pointed out explicitly by Lopez and 
FradkinEll, there are many inequivalent ways to attach 
flux, where particles with like spin see each other as car- 
rying 2ri units of statistical flux, while those with unlike 
spin see each other as carrying 2r 2 units. We will work in 
the picture where the number of flux units seen by other 
particles is independent of the spin, and will deal with 
the case of r\ = r 2 = 1 for simplicity. For this case, the 
statistical gauge field remains a scalar in the spin indices, 
and still satisfies 

V x a(r) = ATtp{r) (52) 

The sequence of canonical transformationsElIll goes 
through exactly as before, and we arrive at the final rep- 
resentation of the electron density projected to the low- 
energy subspace, expressed in second quantized form as 

P(0)= ^2 Pn2n 1 {q)X <Jtnin2 (ci) (53) 
cr,n\n 2 

where we define the spin-cr ME operators 

Xa, ni n 2 (q) = ^{xv^q))* d) aui d a ^ 2 (54) 

k\k 2 

In other words the total density is the sum of the f and 
J. charge densities. One can also define the spin-density 
operators projected to the low-energy subspace 

S z {q) = E ( x u2u 1 (q))*Pn 2 n 1 (q)jdl tVl d (T ,u 2 (55) 
S+(q) = E ( 

{q))*Pn 2ni {q)d\ !Ux d itU2 (56) 

S-(Q)= E( {q)Ypn 2 m{ci)d\ >Vi di, V2 (57) 

kik 2 

The projected spin and charge operators no longer com- 
mute with each other, unlike the full spin and charge 
operators. Instead they form a representation (to lowest 
leading order in q) of the LLL spin-charge algebra dis- 
cussed previously in ref.Ej. For our purposes, it is useful 
to define the spin-flip ME operators 

Y+, ni n 2 {<l)= E (av^i (0)T 4,ki d l<»2 ( 58 ) 

ki &2 

Y-, nin2 (q)= E ( x » 2 vii<l)Td\, Ul dt,u 2 (59) 

k\k 2 

We will be dealing exclusively with fully polarized 
states for simplicity, though the method can be trivially 
extended to partially polarized states as well. Let us arbi- 
trarily designate j as the sign of the majority spin. There 
will be an additional Zccman energy of Ez(ni ~ n^)/2. 
Since the hamiltonian commutes with the ^-component 
of the total spin, this energy can be added to the spin- 
flip excitations at the end. We will therefore ignore the 
Zeeman term when computing the dynamics. As in the 
previous subsection, we commute the hamiltonian with 
Y± and find the hamiltonian matrix in the subspace of 
spin-flip MEs. The final result is 



[H,Y-,m im2 (q)} = (e|(mi) - e T (m 2 ))^-, mim2 (g) 
- (N F i(mi) - N F f(m 2 )) E Y -,nm 2 (q)V x (m 1 m 2 ; nxm) 

nin 2 

where V x is the same exchange matrix element as in Eq. 
([l3|). We also define the spin-dependent HF energies as 

e„(m) = JLj J^ v{q )J2(l-2N Fa (n))\p nm (<l)\ 2 

(61) 

where Np a (m) denotes the spin-dependent Fermi occu- 
pation of the single-particle state m. Note that only 
modes which have (Npi(mi) — N F ^(m2)) ^ need be 
considered. In the fully polarized states Npi(m) = for 
all m, and this restricts us to consider Np^(m) = 1. The 
"negative energy" MEs appearing in the previous subsec- 
tion are absent here, and the TDHF hamiltonian matrix 
is hermitian. 

Let us briefly note that the formalism developed in this 
section can be directly applied to the very interesting case 
of double-layer quantum Hall systems,, which have been 
the subject of much recent interestBEl 

We are now ready to present the results. 



IV. RESULTS 



A. Gaps 



In a previous publication Shankar and the present au- 
thor have presented results for the transport gap in a 
few fractionsty, calculated using the methods of section 
II. We pause here to view the gaps in a slightly different 
light. In many field-theoretic approaches to computing 
the ME dispersions it is assumed that the interactions 
produce CF-LLs whose energy dependence on the CF- 
LL index n is cyclotron-like 

E n = nuj* c (62) 

Indeed, in the CS approach of Lopez and FradkinEl, this 
is the exact result in lowest order, with the bare mass 
m setting the scale for lu*. This seems to be a natural 
assumption. On the other hand, in the LLL, all splittings 
are the result of interactions, and there is no a priori 
reason why they should follow Eq. (|6^) . In our approach 
we can compute these energies, and can therefore test 
this hypothesis. In Table 1 we present the energies of 
the CF-LLs for v = |, for n max — 5. We have set the 
energy of the n = CF-LL to 0, and normalized all 
energies in units of E n= \. The cyclotronic prediction for 
this normalization is E n = n. 
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A 


Ei 


E 2 


E 3 


E 4 


E 5 


0.0 


1.00 


1.75 


2.37 


2.88 


3.32 


1.0 


1.00 


1.48 


1.76 


1.90 


1.89 


2.0 


1.00 


1.32 


1.50 


1.54 


1.43 


3.0 


1.00 


1.22 


1.35 


1.38 


1.27 



It is quite clear that there are strong deviations from 
the cyclotronic form. These deviations become greater as 
the sample thickness increases. In particular, for A > 1.0 
we find the strange result that E$ is less than E4. This 
is another manifestation of the unphysical nature of the 
higher CF-LLs. For higher A there is a distinct tendency 
for the energies to fall towards E\ , to form almost a con- 
tinuum above E\. All these features are present for the 
other fractions as well. 

To emphasize the collective nature of the energies, we 
present the spin-dependent HF energies in the absence of 



Zeeman interactions for v 
in Table II. 



at A = 1.0, with n„ 
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FIG. 1. Spin-polarized ME dispersions for v — 1/3 at 
A = 1.5 for different values of Umax- All energies are in units 
of e 2 1 el. The solid line is the energy of the "naive" ME, which 
corresponds to the strong-field limit in the IQHE. 



n 





1 


2 


3 


4 


E n] /K 2 


-0.36 


0.81 


2.62 


3.39 


3.40 


Enl/K 2 


2.35 


3.74 


4.68 


5.20 


5.07 



Here K 2 is an example of the constant K p defined in 
Section IIIB. The majority (j-spin) and minority (|-spin) 
energies are quite different for the same CF-LL index, 
even when the Zeeman coupling is zero, because the ma- 
jority spin energy has a large negative exchange contribu- 
tion, while the minority spin energy has only a Hartree 
contribution. Neither of the energies has a cyclotronic 
form. 

We therefore come to the conclusion that the cy- 
clotronic form of the energy should be regarded with 
skepticism. 



Let us now turn to the evolution of the dispersions with 
q and n max . As q increases the energy asymptotically 
tends to the transport gap, which itself increases with 
n-max- The gap at n max = 5 is about half its value for 
n m ax — 00. However, the gap at n max = 5 is about 20% 
greater than the true gapEJ for A = 1.5. 
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B. Spin-Polarized Magnetoexcitons 



Figure [l] shows the ME dispersions for v — |, at a 
A = A/Iq = 1.5. We are using a cutoff on the number 
of CF-LL's, n max , denoted by "nx" in the figure. The 
solid line denotes the "naive" ME. In other words, we 
take the state created by the ME operator Xio(q) and 
just compute its energy as a function of q. The energy of 
the naive ME corresponds to the strong-field limit for the 
IQHE. It is clear that including the other ME's, and the 
interaction between them, drastically alters the disper- 
sion. In particular, the magnetoroton minimum present 
in the other dispersions is missing in the naive one. We 
conclude that the strong-field limit is not justified in the 
FQHE. 



0.000 1 ■ 1 ■ 1 ■ 1 ■ 1 

0.0 1.0 2.0 3.0 4.0 

q 

FIG. 2. Spin-polarized ME dispersions for v — 1/3 at 
A = 2.0 for different values of n max - 

The roton minimum is very shallow and the value of 
the energy at q — is quite small. The q at the minimum 
is about I.IIq 1 , slightly smaller than the corresponding 
value for the pure Coulomb interaction. This, is qualita- 
tively consistent with the results of the SMAlj. However, 
for q = the quasiparticle and the quasihole are right on 
top of each other, and our charge density operator may 
not have enough short distance structure to get the right 
numbers for small A. For large A, however, the poten- 
tial becomes very flat at short distances, and we expect 
better results. 

Fig. H shows the same picture for A = 2.0. As can be 
seen, the fractional change in the transport gap as n max 
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increases is even smaller here (the n max = oo value of 
the gap is now about 50% higher than thatrfor n max = 5, 
which is itself 20% higher than the true gapc3). The mag- 
netoroton minimum is even smaller and seems to shift 
slightly towards smaller q. 
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FIG. 3. Spin-polarized ME dispersions for v 
A — 3.0 for different values of max- 



1/3 at 



To continue the trend towards greater thickness we 
show the corresponding results for A = 3.0 in Fig. ||, for 
which the n max — oo value of the gap is is 35% higher 
than that for n max = 5, which is itself within 10% of the 
true gap. The minimum at finite ql has disappeared and 
the minimum is now at q — 0. 
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FIG. 4. Spin-Polarized ME dispersions for v — 1/3 at 
A = 0.0 for different values of n max - Note that the dispersion 
shows an instability for n max — 3, but none for n max > 3. 

Fig. [| shows the same calculation done for the pure 
Coulomb interaction, A = 0.0. For n max = 3 one finds 
that the dispersion becomes unstable for a range of q 
between l.lZjj" 1 to 1.35/q . However, the instability dis- 
appears if n max is increased, and now there is a deep 



magnetoroton 
correct positio: 

ergy are much larger than the true values 



around qlo = 1.25, which is the 
s ' However, the values of the en- 
Thus, while 



a qualitatively correct picture is obtained for n max > 4, 
we cannot trust the numbers for the Coulomb interac- 
tion. The main reason is probably because, as pointed 
out above, our p misses the short-distance physics which 
is important for the Coulomb interaction. 

Since our approach is a dynamical theory of CFs, we 
can obtain all the ME dispersions for a given n rnax . To 
emphasize this, Fig. |^ shows all the ME dispersions for 



A = 2.0, and n ri 
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FIG. 5. All the ME dispersions for v = 1/3 at A = 2.0 for 
n m ax = 4. The lowest mode is well-separated but the others 
form a quasi-continuum. 
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FIG. 6. Fraction of spectral weight carried by the lowest 
mode for v — 1/3 at n max = 4 for different values of A. 
Note that the lowest mode carries vanishing spectral weight 
for large and small q, and that its spectral weight decreases 
drastically with increasing A. 

In Fig. U we show the fraction of total spectral weight 
in the lowest mode for different values of A. The value 
of n max has been set at 4 for all the curves. Two fea- 
tures stand out: Firstly, the lowest mode carries vanish- 
ing spectral weight at small and large q, regardless of A. 
This means that the SMA will overestimate the gap for 
q = and for large q. Secondly, it is clear that for small 
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A, much of the spectral weight is in the lowest mode for 
intermediate values of ql. These twOri-acts are known for 
A = 0.0 from exact diagonalizationsc3. Our calculation 
actually underestimates the weight for A = 0.0 (which, 
to the author's knowledge, is the only case where it has 
been previously computed). However, our results also 
show that the lowest mode carries a decreasing fraction 
of the spectral weight with increasing A, which means 
that the SMA will become less useful for larger A. 
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FIG. 7. Spin-polarized ME dispersions for v = | at A = 1.5 
for different values of n max . 
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FIG. 9. Spin-polarized ME dispersions for v = | at A = 0.0 
for different values of n max . 



Fig. ( J10| ) shows all the dispersions for A = 2.0 and 
n m ax = 8. It is clear that beyond the lowest ME mode, 
there is almost a continuum of excitations (recall that 
these are by no means all the excitations; these are only 
the single quasiparticle-single quasihole excitations). 



The next set of five figures presents results of similar 
calculations for v = | . Here n max = 5 shows instabilities 
for A = 1.5, 2.0, while n max > 6 shows the right behavior 
(we will comment in detail on instabilities in the next 
subsection). The minima appear at qlo = 0.85,2.4. 
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FIG. 8. Spin-polarized ME dispersions for v = | at A = 2.0 
for different values of n m ax- 

While the position of the first minimum is con- 
sistent with the results of exact diagonalizations, the 
second minimum, seems to —he., too high in q (exact 
diagonalizationsEZI and otherEj'LJ numerical results indi- 
cate the second minimum at qlo = 1.6). 
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FIG. 10. All the ME dispersions for v = § at A = 2.0 for 
n m ax = 8. The lowest mode is well-separated but the others 
form a quasi-continuum. 



In Fig. [ll| we present the fraction of the spectral 
weight in the lowest mode for n max = 7. It is clear 
that this is lower for all q at A = 0.0 than in the case of 
i. Once again, as A increases the fraction of the spectral 
weight in the lowest mode decreases. 
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FIG. 11. Fraction of spectral weight carried by the lowest 
mode for v = § at n max = 7 for different values of A. Just as 
in the v = ^ case the lowest mode carries vanishing spectral 
weight for large and small q, and its spectral weight decreases 
with increasing A. 



FIG. 13. Spin-polarized ME dispersions for v 
A = 0.0 for n max = 10. 



at 



Finally, we turn to v — | , for which we present results 
for n m ax — 10 only, and for values of A = 0.0,1.0,2.0. 
We find that the position of first,minimum agrees with 
that from estact diagonalizationstll and CF-wavefunction 
calculations^! (ql rs 0.6), while the second and third min- 
ima seem to be displaced to higher q, just as in the v = | 
case. 
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FIG. 12. Spin-polarized ME dispersions for v = | at 
A = 1.0 and 2.0 for n max — 10. 



Fig. [14] shows the spectral weight in the lowest mode, 
which looks very similar to the | case. We do not re- 
ally trust our results quantitatively for small A. There- 
fore, while one can definitely say that for all fractions the 
spectral weight decreases as thickness increases, it is dif- 
ficult to see any other clear trend (such as a trend with 
increasing denominator) . 
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FIG. 14. Fraction of spectral weight carried by the lowest 
mode for v — | at n max = 10 for different values of A. Just as 
in the v = | case the lowest mode carries vanishing spectral 
weight for large and small q, and its spectral weight decreases 
with increasing A. 

In conclusion, our results with a cutoff on n rnax appear 
to be in qualitative agreement with numerical results for 
A = 0. However, there are quantitative discrepancies 
which make them unreliable for small A. The reasons for 
the failure of our theory to reproduce numbers at small A 
are well-understood, and have been repeatedly mentioned 
above. 

On the other hand, we expect our theory to be quan- 
titatively gaod at large A. Recall that for 2.0 < A < 3.0 
our resultsEll for the transport gap were within 30% of the 
numerical values£J. What needs to be done is to compare 
our results to numerical ones at large A. Unfortunately, 
exact diagonalizations and other numerical methods have 
so far concentrated on the pure Coulomb interaction, and 
few results on ME dispersions are available for large A. 
In the absence of exact results, our results are predictions 
for the ME dispersions. 

We now turn to a brief discussion of instabilities. 
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C. Instabilities Inferred from Magnetoexciton 
Dispersions 



The use of ME dispersions to study instabilities was 
pioneered by Girvin, MacDpaald, and Platzman in their 
seminal paper on the SMAE3. While they did not find 
actual instances where the ME dispersion became nega- 
tive, they did find that for v = 1/7, the magnetoroton 
minimum nearly touched zero. On this basis, they con- 
jectured that in the presence of disorder, the gaps would 
be reduced, and since the potential instability is at finite 
wavevector, that there woijild be a quantum phase, tran- 
sition to a Wigner crystal^. Kamilla and Jainc^l have 
demonstrated within the CF-wavefunction approach that 
for v — 1/9 for the pure Coulomb interaction, the ME 
dispersion indeed becomes negative at a finite wavevec- 
tor, signalling an instability of the incompressible liquid. 
The transition is expected to be first order, and probably 
occurs before (i.e. for v greater than the v for which) the 
ME dispersion becomes negative. 

Our interest will be in instabilities as a function of 
the thickness parameter A. While one can always make 
a CF variable by fiat, by attaching two statistical flux 
quanta to an electron, the experimental reality of the 
CF depends on energetics. It is helpful to think of the 
case_.of v = | in terms of the Haldane pseudopotentials 
I/mE9. These pseudopotentials give the energy of a pair 
of electrons in the LLL with relative angular momentum 
in. The Laughlin state has no pairs in a relative m = 1 
state. The reason this is such a good trial state for the 
Coulomb potential is that V\ is much bigger than V3, and 
there are no other states without pairs in m = 1. This 
is also why the state is incompressible: Every excitation 
must necessarily have an amplitude for creating a pair 
in an m — 1 state, and costs a large energy of order V\. 
Now let us make the sample thicker. As A increases the 
ratio V3/V1 increases, and at some critical value of this 
ratio^he system will make a transition to a compressible 
statecll. This transition is expected to be first ordercH, 
but one may still look for signals of this transition in the 
incompressible state. 

In this spirit, I have analyzed the ME dispersions as a 
function of both A and n max . First consider the v = 1/3 
state. I find that while there appear to be instabilities 
at finite q for small n max (an example is Fig. ^), they 
disappear when one increases n max . So I find no reliable 
signs of instability for v = i. However, for v = |, I 
find that at X cr = 2.65 the ME dispersion becomes zero 
at gi k 0.4. I have increased n max to 10 (much larger 
than the rough cutoff of 2p + 1 = 5 from section), and 
observed that the instability persists. Similarly, I find 



that for v = |, there is a finite wavevector instability at 
A cr = 2.4. It is interesting that the A cr is not the same 
for for fractions which share the same number of flux 
attached, but decreases with increasing denominator. 

Experimentally, such instabilities would be seen as the 
lack of a fractionally quantized plateau beyond a certain 
sample density, even at zero temperature. Presumably, 
the above numbers for A cr are overestimates for the ex- 



perimental values of the thickness parameter, since dis- 
order, which has not been taken into account here, will 
further reduce the gaps. 

Let us now turn to the spin- flip ME dispersions. 



D. Spin-Flip Magnetoexcitons 



Figure [l5| shows our results for v — | for the pure 
Coulomb interaction in the absence of a Zeeman term. 
What is shown is the lowest branch, which corresponds 
in this case to the spin-wave mode. In accordance with 
Goldstone's theorem, the energy vanishes as q — ► 
(quadratically, as it should for a ferromagnet). Upto 
about ql<j w 1 the dispersion for n max — 3 agrees quite 
well with— the exact diagonalization results of Nakajima 
and AokitS, and thus also gives approximately the correct 
spin stiffness. 



0.20 



0.10 



0.00 




nx=3, naive ME 

nx=3 

nx=4 

nx=5 



0.0 



1.0 



2.0 

q 



3.0 



4.0 



FIG. 15. Spin-wave dispersion for v = | at A = 0.0 for 
different values of n max . As Umax increases the separation 
between levels increases, decreasing the interaction between 
naive MEs. 

Once again, the "naive" spin-flip ME is higher in en- 
ergy than the true mode by almost a factor of two in 
this region, and interaction between the various naive 
modes results in closer agreement with exact results. In 
this context the accuracy of the results of Nakajima and 
AokiCj based on flux attachment appear puzzling, since 
the strong-field limit is implicit in their approach. On 
the other hand, A = 0.0 is the worst case for our theory, 
and perhaps it is our agreement with the spin-wave en- 
ergy that is more surprising. It is difficult to compare our 
results to Mandal's worko, since [iftr-the spirit of earlier 
work in the fermionic CS approacriI3'E3 a phenomenologi- 
cal effective mass appears, and the cyclotronic hypothesis 
is used for the energies. 

As q — > 00 the ME energy tends asymptotically to the 
spin-reversed gap (which is considerably overestimated 
by our theory at A = 0.0). It is important to note that 
all higher spin-flip modes are well separated from the 
lowest one for small q, which justifies identifying the low- 
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energy dynamics of the_polarized ^ state with that of the 
quantum ferromagnetcil. 
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FIG. 16. Spin-wave dispersion for v — ^ at A = 1.0 for 
different values of n max - 

We follow the progression as A increases in Fig. [If]. 
Once again, it is unfortunate that numerical results are 
not available for large A, where we expect our theory 
to be quantitatively accurate. Note, that for the spin- 
reversed gao, our prediction is about 30% higher the ex- 
act answerE3 for A as low as 1.5. 



0.10 



0.05 



0.00 



nx=5, lowest mode 

nx=5 next mode 

nx=6, lowest mode 

nx=6, next mode 



2.0 

q 



FIG. 17. Lowest and next lowest spin-flip modes for u = | 
at A = 0.0 for different values of n max . Note that the energy 
of the lowest mode is negative for a substantial range of q, 
and that there is a roton-like minimum. 

Next we present the spin- flip excitations for |, and 
here there is a new feature. While there is always a 
gapless quadratically dispersing mode, it is not always 
the lowest mode. For n ma x — 5 the gapless mode has 
negative energy for small q, while for n max — 6, the gap- 
less mode has positive energy, but the negative energy 
mode persists, and becomes gapped. The lowest mode 
at large q is one which takes an up-spin particle from the 
n = 1 CF-LL and replaces it as a down-spin particle in 
the n = CF-LL. At large q the energy of this excita- 



tion is positive, making the spin-reversed gap positive. 
However, as shown in Figs. [It] and [l^, there is always 
a negative energy mode at small q. In the absence of 
a Zeeman coupling this implies an instability of the pp* 
larized ground state. It is known from previous workcil 
comparing ground state energies that the spin-singlet | 
state is lower in energy than the spin-polarized state for 
small Ez- The zero-temperature transition between the 
states is expected to be first order, but nonetheless, one 
can look for signatures of the instability in the polarized 
state. 
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FIG. 18. Lowest and next lowest spin-flip modes for v — | 
at A = 1.0. All the features observed for A = 0.0 are present. 



In fact, from Fig. [l8| we can infer that at A = 1.0 
(presumably close to the physical range of sample thick- 
ness), the transition to the polarized state should take 
place when a uniform shift of the ME dispersion up- 
wards by Ez makes the dispersion .always positive, i.e., 
for -Mfzr > 0.01. Experimentally^, it is found that 



E 

E z w 0.015e 2 /d. 

In the presence of a Zeeman energy E z capable of sta- 
bilizing a polarized ground state, our result implies that 
the lowest spin-flip mode will have a roton-like minimum 
at finite g, and that the smallest gap to excitations is 
lower than E z . This in turn would indicate a much faster 
reduction of the magnetization with temperature than 
in the case of -|, and a finite-temperature spin-structure 
factor peaked at nonzero momentum. Perhaps the most 
significant consequence of this result is that one can no 
longer think of the polarized | state as a simple quantum 
ferromagnet. 



V. CONCLUSIONS 



I have used the dynamical theory o£-. Composite 
Fermions developed by Shankar and myseliiSEiil to com- 
pute the magnetoexciton dispersions for a fe w gap ped 
fractions. The TDHF method of MacDonalaMla was 
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used to do these computations. Let me reiterate that 
this method is the minimal one in the present case since 
there is only one energy scale. My primary motivation is 
to see how far our theory can take us, and what its limita- 
tions are. To this end, the constraint was imposed in its 
crudest form, by cutting off the CF-LLs at n max . Some 
such cutoff is certainly necessary on physical grounds, as 
explained in section II. 

Our results for spin-polarized ME dispersions apjr-Ur 
qualitative agreement with pH£jt diagonalizationsEKEH, 
and other numerical methodsE3'E3, for the pure Coulomb 
interaction (A = 0.0). This is the worst case for our the- 
ory, since our expression for the charge density is valid 
at long wavelengths, and the Coulomb interaction has a 
significant short-distance part. As a consequence, there 
are significant numerical discrepancies for A = 0.0. 

On the other hand, as the sample thickness increases, 
the interaction becomes softer at high momenta, and our 
theory should become more accurate. To model this, 
we considered a particular interaction with a parameter 
A = \l related to sample thickness. I have presented 
a number of results at larger A. In particular, the cut- 
off dependence becomes negligible at large A, increasing 
our confidence in the results. Generally, we find that 
while the structure of minima found in the ME disper- 
sions for the pure Coulomb case (A = 0.0) persists at 
larger A, the minima become shallower (and even disap- 
pear for ^ around A = 3.0), while the entire ME disper- 
sion moves down in energy. We also find that the frac- 
tion of the spectral weight in the lowest mode decreases 
with sample thickness. Finally, we find that the | in- 
compressible state should become unstable somewhere 
around A cr = 2.65, while the | incompressible state 
should disappear near A cr = 2.4. The instabilities are at 
finite wavevector, and one may speculate that the tran- 
sition might be to a crystal phaseeJ. These are estimates 
for a clean system, and would be revised downwards for 
systems with disorder. Unfortunately, results from ex- 
act diagonalizations or other essentially exact methods 
are not currently available to test our results for thick 
samples. 

Turning to the spin-flip MEs, we find good agreement 
between our results and exact diagonalizationse3 for the 
case v = i, even in the worst case A = 0.0. We find 
that for |, the spin- wave mode is the lowest one, and is 
well separated from the other modes for small q, justifing 
its treatment as a quantum ferromagnctcil. However, for 
v = I we find that the lowest mode is not the spin- 
wave, but a spin-flip excitation that creates a hole in the 
n = 1 CF-LL, and puts a minority spin particle in the 
n = CFpLL. This implies that a quantum fcrromagnet 
treatmentcil that ignores the CFs would be inadequate in 
this case. Furthermore, the lowest excitation has nega- 
tive energy (in the absence of a Zeeman interaction) for a 
substantial range of q, with a minimum at finite wavevec- 
tor. These predictions should be detectable in the spin 
structure factor, and in the temperature dependence of 
the magnetization of the polarized v = | state. 

The most significant deficiency in our treatment is the 
way we handle the constraints. Constraints are known 



to be cruciatJtj'EaOH for compressible systems such 
as v = i. Ideally,^me would like to employ a conserv- 
ing approximation^ in which the constraint would be 
satisfied, such as the one carried out by Read for the 
compressible statescS. An LLL formulation for all the 
principal fractions valid at all q, which is algebraically 
consistent, and which agrees with our approach at small 
q, has recently been postulated by ShankaiO. This for- 
mulation is amenable to a conserving approximation, and 
would be the way to settle matters of principle. However, 
we presented several lines of reasoning to suggest that in 
using our p we are already building in many nontrivial 
renormalizations at tree level, and that the effect of cop-, 
straints would be minimal. The formulation of Shankarta 
shares this philosophy, and provides a partial justification 
of our truncation of p (in the smallness of the higher order 
terms of p pre fcB) . 

In conclusion, our dynamical theory of Composite 
Fermions gives a satisfactory approximate account of the 
low-energy excitations around the gapped FQH states, 
the magnetoexciton modes. In addition to reproduc- 
ing previous results we find some surprising predictions. 
There are many directions in which one can take this ap- 
proach. .Eerhaps the most important are the treatment-of 
disorderEj, and an understanding of the crystal phaseO, 
which I hope to return to in future work. 
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